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Abstract 

An accurate theory describing adiabatic following of the dark, nonabsorbing 
state in the three-level system is developed. An analytical solution for the 
wave function of the particle experiencing Raman excitation is found as an 
expansion in terms of the time varying nonadiabatic perturbation parameter. 
The solution can be presented as a sum of adiabatic and nonadiabatic parts. 
Both are estimated quantitatively. It is shown that the limiting value to which 
the amplitude of the nonadiabatic part tends is equal to the Fourier compo- 
nent of the nonadiabatic perturbation parameter taken at the Rabi frequency 
of the Raman excitation. The time scale of the variation of both parts is 
found. While the adiabatic part of the solution varies slowly and follows the 
change of the nonadiabatic perturbation parameter, the nonadiabatic part 
appears almost instantly, revealing a jumpwise transition between the dark 
and bright states. This jump happens when the nonadiabatic perturbation 
parameter takes its maximum value. 

I. INTRODUCTION 

Stimulated Raman adiabatic passage (STIRAP) resulting in the population transfer be- 
tween the states which are not directly coupled [1], electromagnetically induced transparency 
(EIT) via adiabatic following of the dark, non-absorbing state [2], nonresonant pulse exci- 
tation of the two-level atom [3] are just a short list of phenomena in quantum optics where 
adiabatic processes are considered. Generally, one can find in any part of physics problems 
concerned with adiabaticity, which are treated almost similarly. Among them we can men- 
tion multiphoton resonances induced in atoms and molecules by a strong low-frequency field 
[4]- [7], wave-packet dynamics in physics and chemistry or, so-called, "femto-chemistry" [8]- 
[10] and slow atomic and molecular collisions [11]- [22]. 

A fully adiabatic process takes place if the Hamiltonian of the problem is diagonal. Then 
any time dependence of this diagonal Hamiltonian does not change the state of the quantum 
system except its phase. In general, the Hamiltonian, describing the interacting systems, has 
both diagonal and nondiagonal components. In most cases dealing with adiabatic excitation, 
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the problem can be reduced to the consideration of the two-level system perturbed by the 
diagonal and nondiagonal interactions. Then, both components of the Hamiltonian, i.e., 
the two- level splitting and the nondiagonal perturbation, are time dependent. One can find 
instantaneous eigenstates of such a Hamiltonian. These states are called adiabatic states if 
the system, being initially in one of them, follows this state. In contrast, the initial states, 
where the transition takes place, are called diabatic states. 

The diagonalization of the instantaneous Hamiltonian takes into account the diabatic 
states coupling, incorporating the interaction parameter into the instantaneous eigenvalues. 
However, a new coupling appears again because of the time dependence of the transformation 
from the initial, diabatic basis to the adiabatic basis (see, for example, [10]). Inasmuch as the 
coupling between instantaneous eigenstates does not vanish, the following of the adiabatic 
states can never be perfect. Since it is almost impossible to get rid of the state coupling by 
the time dependent transformation of the basis functions, the choice of the adiabatic basis is 
aimed to minimize the transition probability between the states, almost approaching in this 
sense to the fully adiabatic Hamiltonian (i.e., the diagonal one) when the coupling can be 
neglected. The estimation of the fraction of the system that is transferred from one adiabatic 
state to another during the interaction (nonadiabatic correction) is a main problem in this 
approach. This nonadiabatic correction is a measure of the nonadiabaticity of the process. 

It should be obvious that the searching for the adiabatic basis is an approximate ap- 
proach. It was developed for the first time to treat the behavior of the two-level system sub- 
ject to an excitation that couples the states if both the two-level splitting and the excitation 
are time dependent. Only a few models of such a two-level system have analytical solutions. 
Among them are the Rosen-Zener model with H d i ag = const and H nond = Vq sec h(rt), 
the Demkov-Kunike model with H diag = E + E tanh(rt) and H nond = V sec h(rt) or 
Hnond = const, the Landau-Zener model with H diag = \t and H nond = const, and the 
parabolic model with H diag = at 2 + b and H nond = const. Although the last model has not 
an exact solution, it has been studied comprehensively (see, for example, Ref. [23]). Here, 
H d iag is the amplitude of the diagonal component of the Hamiltonian (i.e., the two-level 
splitting) and H nond is the amplitude of the nondiagonal component (i.e., the coupling per- 
turbation). The symbols V , E, E Q , r can represent arbitrary constants and t is time. One 
can find these models and their solutions in, for example, [10]. 

The advantage of these models comes from the possibility to reduce their Schrodinger 
equations to second order differential equations (with time dependent coefficients) whose so- 
lutions are known. For example, the Schrodinger equation for the Rosen-Zener model can be 
reduced to an equation whose solution is the hypergeometric function and the Schrodinger 
equation for the Landau-Zener model can be solved in terms of parabolic cylinder functions. 
For an arbitrary Hamiltonian, the quasi-classical approach based on the Wentzel-Kramers- 
Brillouin (WKB) approximation was developed by Dykhne [24]- [25]. In this method, the 
amplitude of the nonadiabatic transition is calculated under the assumption that the adia- 
batic levels cross in the complex time plane and the main contribution of the nonadiabatic 
coupling is given at the time of the level crossing. The latter is usually the branch point 
of the phase integral of the energies. The result becomes universal, i.e., independent of the 
nonadiabatic coupling model and it is sufficient to take the leading term of the coupling at 
the crossing point. 

In this paper we consider adiabatic following of the dark, nonabsorbing state in the 
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three-level system excited by two resonant fields. The adiabatic following of the dark state 
results in STIRAP (population transfer). It was proposed in Ref. [26] where the so called 
counterintuitive Raman pulse sequence emerged as a result of the search for the generaliza- 
tion of the Liouville equations for the iV-level system using SU (N) coherence vector theory 
[27]- [26]. This pulse sequence consists of two fields coupling two low energy levels 1 and 
2 with one common excited state 3. If one of the levels (for example, 2) is initially empty, 
then applying first the field which couples this state with level 3, and then the field coupling 
the populated state (for example, 1) with 3, it is possible to transfer the population of state 
1 to state 2, without appreciably populating the intermediate state 3. 

Later, the importance of the adiabaticity in the counterintuitive pulse sequence devel- 
opment was realized in [29]- [31]. There is a particular superposition of states 1 and 2 that 
does not interact with the coupling fields and, if the development of the field amplitudes 
in time is properly chosen, this superposition state changes from state 1 to state 2. If the 
three-level atom follows this superposition state, the atom population is transferred from 
state 1 to state 2 by the pulse sequence without populating the intermediate excited state 3. 
This noncoupled superposition state was first introduced by Arimondo [32] to explain qual- 
itatively the dark resonance as population trapping in this state [33]- [35]. The noncoupled 
state is often referred to as a dark state. 

It is obvious that, if the dark state changes in time, a process must exist which tends 
to empty this state. The condition minimizing the dark state depopulation is formulated 
in [29]. This is done ad hoc, without an estimation of the nonadiabatic correction for the 
excited probability amplitude. However, numerical calculations show that, if this condition 
is satisfied, the adiabatic population transfer 1 — > 2 is almost perfect. Some attempts were 
undertaken to find a rigorous justification of the intuitively found adiabaticity condition and 
to estimate the amplitude of the excited state during the STIRAP pulse sequence (see, for 
example, Refs. [36]- [37]). In [36], the so called ramp pulses were considered, which allow an 
exact solution. Furthermore, with the help of the method developed for the two-level system 
by Dykhne (see the discussion above and Refs. [38]- [39]), the nonadiabatic amplitude of 
the excited state 3 is estimated for Gaussian and hyperbolic secant pulses. This is possible 
because the Liouville equation for the two-level atom in terms of the SU (2) coherence vector 
(Bloch- vector model, see, for example, [40]) coincides with the Schrodinger equation for the 
state probability amplitudes of the three-level atom excited by two resonant fields (see, for 
example, [31], [36], [2]). However, the authors of [36] admit that the analytical approxi- 
mations for the nonadiabatic corrections "have been introduced ad hoc without derivation" 
and they "would really like to see more detailed investigations of the analytic behavior of the 
system discussed 1 . 

Fleischhauer with coauthors [37] developed a different approach introducing higher-order 
trapping states. They define a nth-order generalized adiabatic basis, which is similar to 
the superadiabatic basis introduced for the two-level system in [41]- [44]. By successive 
transformations from the diabatic basis to the adiabatic basis, then from this adiabatic 
basis, which can be considered as the adiabatic basis of the first order, to the next, i.e., 
the second order adiabatic basis, etc., the solution is presented as an infinite product of 
transformations. The general expression for the n-th transformation matrix is presented. 
However, it is hard to implement this scheme for an arbitrary pulse sequence. 

In this paper we develop a new method that allows us to calculate the adiabatic and nona- 
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diabatic components of the solution with controllable accuracy. The adiabatic component 
describes the part of the atomic probability amplitude visiting the excited atomic state dur- 
ing the pulse train and coming back to the dark state when the pulses are gone. The atom 
evolution following the adiabatic component of the solution resembles the excitation-de- 
excitation process induced by a soliton in the two-level atom. The nonadiabatic component 
is that part which is lost from the dark state and describes the fraction of the atomic prob- 
ability amplitude that is left excited after the pulse train. We compare our result with the 
previous theories reported in [36] and [37]. We found corrections to the theory presented 
in [36] and show that our result corresponds to the calculation of the infinite number of 
transformations proposed in [37]. As an example of the validity of the method, we present 
also the solution of the Rosen-Zener and Landau-Zener models by our method (see the 
Appendix and the end of Sections VI and VII). 

The paper is organized as follows. In Sec. II we present the general formalism employed 
in the description of the three-level atom excited by two resonant fields. The transformation 
to the basis of the bright and dark states is derived. It is shown that the system evolves 
between two states, i.e., "bright" and "common" (they are specified in Section II). In Sec. 
Ill we consider the time evolution of the atomic state vector for the A-scheme of excitation 
if the field amplitudes are time dependent. Bloch-like equations are derived. In Sec. IV the 
adiabatic following approximation is presented. The adiabatic solution for the stimulated 
Raman adiabatic passage (STIRAP) is found in Sec. V. Nonadiabatic corrections are de- 
scribed in Sec. VI. The case when the Raman Rabi frequency changes in time is considered 
in Sec. VII-IX. 

II. THREE-LEVEL ATOM INTERACTING WITH TWO RESONANT FIELDS. 

GENERAL FORMALISM 

We consider a three-level atom shown in Fig. 1 (a). The arrows indicate the transitions 
induced by the coherent fields. This excitation scheme is classified as the A scheme. We 
define the level that is common for both transitions as 3. The others are designated by the 
numbers 1 and 2, level 2 being of higher energy than 1 and initially not populated. The 
dynamic part of the Hamiltonian of this three-level atom, excited by two resonant fields 
Ei(t) = Ei cos(fiit + <fi) and E 2 (i) = E 2 cos(Q 2 t + ^2), is 

3 

H = ^nPnn ~ (b^P^^ + B 2 P 23 e^ t+ ^ + /l.C.) , (1) 

71=1 

where u) n is the energy of the state n. Planck's constant is set equal one (h — 1) for 
simplicity. The operators P mn are defined by P mn = \m) (n\, where (n\ and \m) are bra and 
ket vectors of the states n and m in the Schrodinger representation. The interaction constant 
(Rabi frequency) B n = (d n3 • E n ) /2 = (d 3 „ • E n ) /2 depends on the dipole-transition matrix 
element, taken real so, d n3 = d 3n , and on the field amplitude E n . The rotating wave 
approximation is taken into account. 

If the fields Ei(t) and E 2 (t) are in exact resonance with the relevant transitions, the 
Hamiltonian, Eq. (1), can be made slowly varying by transforming it into the interaction 
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representation (IR). This representation is defined by a canonical transformation by means 
of the unitary operator 



T = exp i UnPnnt 



(2) 



n=l 



The wave function of the atom in the IR is defined by | $>(£)) —T\^f(t)), where \^{t)) is the 
wave function in the Schrodinger representation. |3>(t)) satisfies the Schrodinger equation 
with the effective Hamiltonian [45] 



H = THT- 1 + iTT-\ 
This Hamiltonian has the explicit form 

H = -5iPi 3 e^ - B 2 P 23 e iip2 + h.c, 

or in matrix notation 



(3) 



(4) 



n 



fiie^ 1 

P 2 e^ 2 

B x e'^ B 2 e-^ 2 



(5) 



The P/wAr-operators are defined for the vectors \M) and (N\ of the interaction representation 
differing from the states \m) and (n\ by the phase factors exp(— iu m t) and exp(iu> n t). 

Assume that the field amplitudes B ± , B 2 and the phases ip±, ip 2 are constant. Then the 
Hamiltonian, Eqs. (4)-(5), is diagonalized by a unitary transformation 



H = QHQ- 1 = yjB* + Bl (P33 - Pit) , 



(6) 



where 



Q 



1 

75 



COS OL 



cos a -v^e^ 1 sin 



a 



sin a 



-up 2 



cos a 



(7) 



and tana = B\jB 2 . The new ket vectors \n) are related to the former ones \N), defined in 
the interaction representation, as follows |$) = Q |<E>). The explicit relations are 



\/2 

|l> = ^— (e^sina |1) + e^ 2 cosa |2) + |3)) , 



(8) 



1 2) = e^ V2 cos a |1) - e-* 1 sin a |2) , 



(9) 



1 3) = y~ (e^ 1 sin a |1) + e llp2 cos a |2) - |3)) . 



(10) 
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The basis \n), in which the Hamiltonian is diagonal, is called the basis of the quasi-energy 
states [46], [47]. It coincides with the basis of the dressed states if the limit of the classical 
field is taken for them. 

States |l) and |3) are mixtures of all unperturbed states |1), |2) and |3), whereas state 
1 2) does not contain the common state |3). If the atom is in state |2), a mixture of the 
ground state sublevels 1 and 2, the atom does not leave this state since it is an eigenstate of 
the interaction Hamiltonian, Eq. (6). The bichromatic field E(t) = Ei(t) + E2(t) does not 
interact with such an atom, and, consequently, it is not excited. Therefore, we call this state 
a dark state and designate it \d). This state was introduced for the first time by Arimondo 
in Ref. [32], who called this state a non-coupled state. He introduced the coupled state as 
well, which interacts with the bichromatic field E(t). Following Arimondo we define the 
state 



|6> 



e lLpl sin a 



|1) + e iV2 cosa |2), 



(11) 



which is orthogonal to the dark state. We call this state a bright state, since for the A- 
scheme, if the atom is in this state, it is excited by the bichromatic field and then the 
luminescence from state |3) may follow. Excitation can take place because the bright state 
is not an eigenstate of the interaction Hamiltonian, Eq. (6). 

The states \d), \b) and |3) are mutually orthogonal and can be chosen as a new basis. 
We call this basis dbc, designating the state |3) by the letter |c), since it is common for the 
induced transitions. Because we will often refer to these states, they are presented below in 
a common set of equations to simplify further citation 



\d) = e~ ilp2 cos a |1) - e-* 1 sin a |2) , 



(12) 



1 6) = e ilpl sin a |l) + e^ 2 cosa |2) , 



(13) 



|c) = |3>. 

The interaction Hamiltonian, Eqs. (4)-(5), is transformed in this basis as follows 

7~Ldbc = SHS 1 = — B ^Pbc + P c bj , 



(14) 



(15) 



or in a matrix form 



T~tdbc — 




00 s 
B 



(16) 



where B = B\ + B\ and the unitary operator of the canonical transformation is 



e llp2 cos a —e lipi sin a 
e -»vi s i n a e ~ lip2 cos a 
1 



(17) 



We call this basis the ci&c-representation. 
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Figure 1 (b) shows the excitation scheme in the <i&c-basis where only the b and c states are 
coupled. The b and c states do not correspond to a defined energy. To present schematically 
the dbc states and the transitions between them, we choose the state population before the 
excitation for their relative position in the diagram. So, the vertical scale in Fig. 1 (b) 
represents the initial population of the states counted from the bottom to the top. 

The Hamiltonian in the <i6c-basis, Eq. (15), resembles the interaction Hamiltonian of 
the two-level system be, excited by one resonant field with an effective interaction constant 
B. The dynamic evolution of the atom is described by the Schrodinger equation with this 
Hamiltonian containing only the transition operators between states b and c. If, for example, 
initially the d and b states are populated and c is not, 

\&bc(0)) = C d \d)+C b \b), (18) 

the dynamics of the three-level atom follows the very simple formula 

\$ bc (t,±)) = C b 

where xr — 2-B is the Raman Rabi frequency. This Rabi frequency is defined as twice the 
interaction constant in the Hamiltonian, Eq. (15), according to the conventional definition 
adopted for the two-level system [40]. The initial condition specified in Eq. (18) has a simple 
relation to the IR basis 

|$ 6c (0)) = (C b e ilpl sin a + C d e~ iip2 cos a) |1) + (C b e ilp2 cos a - C d e~ iipi sin a) |2) . (20) 

If the initial condition is |<& 6c (0)) = |1), then C b = e~ tlfl sin a and C d = e llf2 cosa, which 
means that the phase relation between the amplitudes of b and d states is fixed as C b = 

C de -*(<pi-<w)tana. 

Suppose that the pulse duration of the bichromatic field, t2 n , is chosen such that the 
effective two-level system be experiences a so called 27r-pulse, i.e., Xr^-k = 27r. During the 
first half of the pulse, the population of level b is completely transferred to state c. During 
the second half of the pulse, the population of the excited state c is transferred back to state 
b. The state vector of the atom at the end of the pulse is 

\®bc(t2n)) = -C b \b)+C d \d). (21) 

Equation (21) shows that the 27r-pulse results in a n phase shift of the b state probability 
amplitude. Usually, for the simple two-level atom, this ir phase shift of the ground state 
after a 27r-pulse is not revealed in the observables. In our case, apart from the two levels 
b and c, we have the third level d whose phase is well defined with respect to the phase of 
state b. The phase definition is determined by the phase relation of the components of the 
bichromatic field [see comments after Eq. (20)]. In section IV we show the importance of 
this point. 

III. TIME DEPENDENT AMPLITUDES 

If during the pulse excitation the field amplitudes E 1; E 2 are time dependent and they 
have the same time evolution satisfying the condition Bi (t) / B 2 (t) = tana = const, then 



cos 



XRt 



\b) + i sin 



XRt 



+ C d \d), 



(19) 
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the developments of the \d), \b) states in the ket-vectors |1), |2) have constant coefficients 
and the S'-transformation is time independent. In this case the solution of the Schrodinger 
equation is simple and the development coefficients of the state vector \Q dbc ), 

\$ dbc (t)) = C d \d) + C b \b)+C c \c) : (22) 

are C d (t) = C d (0), C b (t) = C b (0) cos[0(f)/2] and C c (t) = iC 6 (0) sin[0(t)/2], where C d (0), 
C b (0) are the initial values of the probability amplitudes and 0(t) = 2 y B{r)dT is the 
pulse area of the bichromatic field. This case corresponds to so called matched pulses [48]- 
[49]. 

If there is a time shift T between the pulses Ei(t) and E2(t), one can again introduce the 
bright and dark states, employing the time dependent ^-transformation. Then the evolution 
of the state vector of the atom in the <i&c-basis is given by 

\* dbc {t)) = S{t)\Q{t)), (23) 

where | is the state vector in the interaction representation. Taking the time derivative 
of the equation (23), one can obtain the Schrodinger equation 

^^ = -<W*c|$« Bc > (24) 

with the modified Hamiltonian 

7~Ld bc = 7~Ld bc + i S S~ x , (25) 
where Tid bc is defined in Eq. (15) and 



i S S' 1 = i a 



(26) 



The first part, Tt d bc, of the modified Hamiltonian induces transitions between the \b) and 
|c) states with the rate xr — 2-B, and the second part, % S induces transitions between 
states \d) and \b) with the rate 2 a (see Fig. 2). 

The development coefficients of the state vector \& dbc ) satisfy the equations 

Z d = -aY b (27) 
Y b = -BX C + a Z d (28) 

X c = BY b (29) 

where the substitution Z d = C d exp(— i(p2), Y b = C b exp(iipi), X c = —iC c exp(iipi) is made to 
deal with real numbers. Here the phases tpi and v?2 are assumed to be constant throughout 
the excitation process. Equations (27)-(29) remarkably coincide with the Bloch equations 
for an abstract two-level system g — e excited by a field with frequency cu slightly tuned from 
resonance (g and e being ground and excited states, respectively) [2]. Expressed in terms 
of the Bloch-vector components, which are the following combinations of the density matrix 
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of the two-level system p: u — iv = 2p ge exp(—iujt) and w = p gg — p ee , these equations are, 
Ref. [40], 



w= -xv, 



(30) 



v— —Au + x w . 



(31) 



u= Av. 



(32) 



Here A is the detuning from resonance and x is t ne Rabi frequency. If one makes the 
substitution w = Zj, v = Y&, u = X c for the variables and x A = 5 for the parameters, 
both sets of equations are identical. 

The coincidence of the equations enables us to use the adiabatic following approach de- 
veloped by Crisp [3] to describe nonresonant excitation of the two-level system. A similar 
approach was applied by Laine and Stenholm (LS) [36] based on the ideas of the adiabatic 
following developed by Dykhne [7], [25] for the two-level system with time dependent split- 
ting and coupling parameters. In the LS approach, [36], the instantaneous eigenstates of 
the three-level atom excited by two resonant pulses are found. This basis is called the adi- 
abatic representation. The transformation to the instantaneous basis [see Eq. (7)] and the 
instantaneous Hamiltonian diagonal in this basis [see Eq. (6)] are time dependent. There- 
fore, the adiabatic states coupling, i Q Q~ l , also appears in their consideration. In spite of 
being different in structure, our Schrodinger equation in the changing rffec-basis and the LS 
equation in the instantaneous eigenstate basis can both be reduced to the equation for the 
two-level system. However, our equations (27)-(29) are in a one to one corespondence to the 
Bloch equations while the LS equations match the equations for the two-level density matrix 
elements p eg , p ge and p gg — p ee . Although this difference is not crucial, it brings, however, 
some convenience in our case, because our equations are expressed for real quantities. 



Crisp gave in his paper [3] a qualitative and quantitative description of the adiabatic fol- 
lowing approximation for the two-level atom excited by a nonresonant field. The qualitative 
explanation employs the Bloch-vector model as follows. The Bloch equations can be written 
in the form 



if the relaxation times 7\ and T 2 are very long compared with the pulse duration. The Bloch 
vector R is given by 



IV. ADIABATIC FOLLOWING APPROXIMATION 



R= H e // x R, 



(33) 



R = Xe x + Ye y + Ze 



(34) 



and the effective field H e // is 



H e// = -X^x + Ae 2 . 



(35) 
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Equation (33) has the geometric interpretation that the Bloch vector tries to precess about 
the effective field as H e // varies both in magnitude and direction. If the condition |A| ^> \x\ 
is satisfied, the precession frequency of the Bloch-vector, a/x 2 + A 2 , will be large compared 
with the rate of change of the field x(0- ^ initially the Bloch-vector points parallel or 
antiparallel to e z [X(0) = Y(0) = and Z(0) = ±1], then precessing with high frequency 
|A| about the effective field, the Bloch-vector remains nearly parallel or antiparallel to the 
effective field as it moves adiabatically. Qualitatively, adiabatic movement of the vector 
tt e ff means that the Bloch-vector does not slip off from the direction of the effective field 
if the rotation of the Bloch-vector around the H e //-field is much faster than the change 
of the H e /_/— field direction. Quantitatively, this sets a condition that the precession period 
T p ~ 2tt/ I A | is to be much shorter than the time scale of the effective field variation. In 
this limit, the angle between the Bloch vector and the effective field will remain very small. 
The deviation of the Bloch vector from the effective field is described by the nonadiabatic 
corrections. In the next section we consider how small this deviation angle is for particular 
pulses. 

The approximation of the adiabatic following allows also a quantum mechanical inter- 
pretation. The Hamiltonian in the dbc time varying representation, Eq. (25), consists of two 
parts. One part couples two, initially unpopulated states, b and c. The coupling parameter 
B is time varying. The other part couples the populated state d with the unpopulated state 
b. The coupling parameter a is also time varying. In order that the atom is subject to the 
adiabatic following of the changing d state (the atom stays in this state), one has to choose 
the pulse train with the Raman Rabi frequency B much larger than a when the latter takes 
its maximum value [2]. Also the time variation scale of a must be much longer than the 
Raman Rabi period. These conditions can be made plausible by the following speculations. 

Since initially only state d is populated, the atom starts to evolve when a (t) becomes 
nonzero. When the slow process a brings a small fraction of the population to state b, the 
fast Rabi oscillation via the 27r-cycle (see the discussion at the end of Section II) changes the 
sign of the lb-amplitude from plus to minus, reversing back the population transfer process 
from d — > b to b — > d. As a result, the population of the 6-state does not increase and 
remains small. One can roughly estimate its mean value by solving the equations (27)-(29) 
with the initial condition X c (0) = Yb(t) = 0, Zd(0) = 1 and assuming that xr an d ot are 
constant, i.e., neglecting the derivatives of these parameters, 

B 2 a 2 I o 

Z d{t) = 72 + 72 C0S V B2 + a t i ( 36 ) 

B 2 + a B 2 + a 

Ct I 9 

Y b (t) = . sin V B 2 + a t, (37) 

VB 2 + a 

X c (t) = ^ B 2 ^1 - cos \Jb 2 + a?j . (38) 

The solution shows that the change of the amplitude of the d-state is small if B = 
a/ B\ + B 2 3>a. This is exactly the condition of the adiabatic following formulated in 
Ref. [29] for STIRAP. 
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It should be emphasized that the higher time derivatives of the mixing parameter a(t), 
neglected above, play a crucial role in the adiabatic population transfer. However, the 
solution (36)-(38) gives some qualitative hint on the choice of the relation between B and 
a. 



V. ADIABATIC SOLUTION 



The adiabatic following approximation can be applied to the consideration of STIRAP 
and electromagnetically induced transparency (EIT) because in both cases the atom follows 
the dark, noncoupled state. A first attempt to study adiabatic following in EIT was under- 
taken in Ref. [2]. In this section we consider the application of this approach to STIRAP. 

The stimulated Raman adiabatic passage assumes that before the application of the Ei(t) 
and E 2 (t) pulses the atom is in state |1). The duration of the excitation as well as the pulse 
sequence must be chosen such that at the end of the pulse train the atom is left in state |2). 
It is expected that during this process the atom stays in the dark state. However, this state 
itself changes since the coefficients cos a and sin a of the development of the dark state in the 
vectors |1) and |2) [see Eq. (12)] change in time. The parameter a rises from zero to ir/2, 
so \d) = |1) before the excitation and \d) = —e~ l ^~vi) |2) after it. Since tana = B 1 /B 2 , 
the condition imposed on a means that the _Bi(t)-pulse must be delayed with respect to the 
B 2 (t)-pu\se. We have to emphasize that the phase of the final state 2 after the pulse train 
is related to the phase of the initial state 1 according to Eq. (12). 

If we choose two identical, bell-shaped, delayed pulses having a hyperbolic secant shape, 
then the interaction constants evolve in time as follows 



B n (t) = B sech[r(t-t n )], (39) 

where n = 1 or 2, t n is the time when the n-pulse has maximum amplitude, and r is the rise 
and fall rate of the pulse edges. The mixing parameter a(t) increases monotonously if the 
condition t\ > t 2 is satisfied. This pulse sequence was considered by Laine and Stenholm in 
Ref. [36]. 

Let us analyze the constraints imposed on the parameters r, ti and t 2 to have a changed 
from zero to ir/2. The mixing angle a varies according to the relation tana(i) = Bi(t) / B 2 {t) . 
For the pulse sequence specified above, an explicit form of this relation is 

l + Dtanh[r(f-f )] , , 

tanat = -f-7 ^7, 40 

w l-Dtanh[r(t-f )] 

where D = tanh(rT/2), t = {t\ +t 2 ) /2 is the mean time of the maxima and T — t\ — t 2 
is the time interval between the maxima of the pulses. Suppose that, initially, the atom 
is in the ground state 1 and we start the atom evolution from the initial mixing pa- 
rameter ««„(— 00) satisfying the condition tan(aj n ) = 0.01, which means that essentially 
Cd{— oo)e~ l</32 = cos(o;j„) ~ 1 and Cb(— oo)e* Vl = sin(o;j ra ) = 10~ 2 . We stop the atom evolu- 
tion at tan[o!/j n (+oo)] = 100. So, if the final state coincides with state 2 then Cd(+oo)e~ tlpl = 
sin(a/ in ) w 1 and Cb(— oo)e^ 2 = — cos(o;/j n ) = — 10" 2 , where the phase change of state 2 
is taken into account. From Eq. (40) it follows that tana(±oo) = (1 ± D)/(l =f D) and at 
the condition imposed on a in and we have rT = 4.6. Since state b is strongly coupled 
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with state c [the coupling is B(t)}, we have to keep the initial population of state b as small 
as possible [C&(— oo) must be close to zero]. Otherwise, the probability amplitude C b (— oo), 
if not infinitely small, will be spread among the b and c states by the pulse train and pop- 
ulation transfer 1 — > 2 will be imperfect. Therefore, on the one hand, the initial value of 
the mixing parameter a must be as small as possible to have complete population transfer 
1 — > 2. For a small initial mixing angle ai n , the relation between oti n and T becomes simple, 
i.e., tanh(rT/2) rjI - 2«j n . The smaller the initial value of the mixing angle a, the larger 
the product rT or the pulse separation T is. On the other hand, if the distance between the 
pulses is large, the value of the coupling B(t) at t becomes small: the larger the distance, 
the smaller the coupling. However, the adiabatic following demands a large coupling B 
at to- Therefore, one has to choose the optimum value of the pulse spacing satisfying two 
conditions simultaneously. The distance between pulses is to be as large as possible to have 
a small value of a in and, at the same time, B(t ) should be kept as large as possible. Below 
we give some arguments how to find this optimum value. 

If rT — > 0, then according to Eq. (40) we have a — > n/4 throughout the excitation 
and the dark state does not change in time. It has the probability amplitude \Cd\ = n/2/2. 
State b, having the initial population |C&(— oo)| 2 = 1/2, is depopulated with the rate B 
[C b (t) = C 6 (-oo)cos(0(t)/2), C c (t) = iC b (-oo)sm(9(t)/2), see the definition of 6(t) in 
Section III]. For small values of rT, the time interval of the a change is small. For example, 
if rT = 0.5, then cos(a in ) = 0.855, sm(a in ) = 0.519 and cos(af in ) = 0.519, sin(aj in ) = 0.855, 
which corresponds to the change of a from cti n = 0.1747T to ai/j n = 0.3267T during the pulse 
train. In this case only the fraction (0.855) 2 = 0.731 of the atomic population is transferred 
to state |2) if the atom adiabatically follows the dark state. Another fraction (0.515) 2 = 0.269 
of the atomic population participates in the process of excitation and de-excitation between 
states b and c. This means that the population transfer via the change of the amplitude of 
the dark state components is possible only if the time interval between the pulses exceeds a 
certain value. We choose the value rT = 5 since in this case the initial population of state 
b is 4.5 • 10 -5 so that we neglect this population in our further consideration. 

Figure 3 (a) shows the pulse train with rT = 5 along with the change of the mixing 
parameter a during the excitation. On the same plot the dependence of the derivative 

rD sec h 2 [r(t — t )] rsinh(rT) 
a ~ 1 + D 2 tanh 2 [r (t - t )] ~ 1 + cosh(rT) cosh[2r(t - t )] 

is shown. This derivative takes a maximum value « max (^o) — r D at t = t . For example, if 
rT = 5, then a max {to) = 0.987r. In Fig. 3 (b) we compare the time evolution of the Rabi 
frequency Xi?/2 = B(t) = y/ B\{t) + B 2 (t) with the evolution of the derivative a (t). 

The Rabi frequency determines the transition rate between states b and c, whereas the 
derivative of the state mixing angle a specifies the transition rate between states d and b 
[see Eqs. (27)-(29)]. Since at t — to the parameter a (t) takes its maximum value and 
B(t) has its minimum, the adiabatic following condition is B(to) (t ) or, explicitly, 
\p2Bo r sinh (rT/2). Figure 4 (solid lines) shows the time dependence of the amplitudes 
X c (plot a), Y b (plot b) and Z d (plot c) obtained by the numerical solution of the equations 
(27)-(29) for the pulse train with t x = 2.5/r, t = 0, t 2 = -2.5/r, T = 5/r, B = 42.8r. 
The relation between B Q and r corresponds to the ratio B(t )/ a (t) = 10. We see that the 
deviation of the Z rf -amplitude from its initial value (~ 0.005) and the maximum absolute 
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value of the Y^-amplitude (~ 0.01) are different from the values predicted by the simplified 
solution, Eqs. (36)-(38). Only the maximum value of X c (~ 0.1) fits the value of the 
simplified solution. 

To estimate the actual values of the probability amplitudes of the dbc states during 
the excitation and find the borders within which the adiabatic following of the dark state 
takes place, we follow the theory developed by Crisp [3] for the case when the condition 
B(to) ^>a (t) is well satisfied. 

In this case we can use the expansion in a power series of the parameter a (t) for the 
solution of the equations (27)- (29). Then the first term of the expansion is found by setting 
Z d (t) = 1 in the equations (28), (29) and then solving them. The solution is 



Y, 



(t) = I cos I Btnjdn 

J— OO J T 



a (r)dr, 



(42) 



X c (t) 



sin 



S(n)dn 



a (t)cZt. 



(43) 



This is the part of the general solution that is linear in a (£). Z d satisfies the equation 

Z d (t) = l- [ a(r)X b (r)dr, 



(44) 



which is the formal solution of Eq.(27). Therefore the change of the dark state amplitude Z d 
is nonlinear in a (t) and can be presented in products of as. The first contribution of a (t) 
to Zd can be found as the square of a, i.e., in the second term of the expansion. Substituting 
the corrected Z& into Eqs. (28), (29) (instead of Zd = 1), one can find the next term in the 
expansion of X b and Y c . Then the substitution of the thus found X b into Eq. (44) gives the 
next term of the expansion of Z d , etc. In this paper we consider only the linear corrections 
to Y b and X c . 

Our solution, given by the equations (42)-(44), differs from Crisp's solution in the B- 
time dependence. He considered the case when B is constant and a is time dependent. Fig. 
4 (a-c) shows the comparison of the numerical solution of the equations (27)-(29) with Eqs. 
(42)-(44). They are indistinguishable and shown by the same solid lines. 

We find the adiabatic and nonadiabatic components of the analytical solution by applying 
two different procedures. The adiabatic components of the equations (42)-(43) are calculated 
by integrating them by parts. For example, the first step of the X c (t) calculation is 



X c (t) 



B{r) 



d ( cos 



Birjdn 



a 
B 



cos 



Birjdn 



dr. (45) 



where (a /B^j T is the derivative with respect to r. Repeating these steps several times, we 
obtain 



Y b (t) = a 2 — «4 + «6 — «8 + 



(46) 



X c {t) = OL\ — 0:3 + 0:5 — «7 + 



(47) 
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where a n =a n -i /B, ao = a and it is assumed that a (t) is a bell-shaped function with first 
and higher derivatives equal zero at t — — oo. To find a similar expansion of equation (44), 
we introduce the function Q(t) = f Q B{r)dr. Then the equation for Z d (t) takes the form 



with 



Z d (t) = l-A c (t)-A s (t), (48) 
A c (t) — / dri a (n) cosQ(ri) / dr 2 a (r 2 ) cosf2(r 2 ), (49) 

J —oo J —oo 

/t /-Tl 
dri a (ri) sinfi(ri) / g?t 2 a (r 2 ) sin fi(r 2 ). (50) 
-oo J —oo 

The functions A c (i) and A,(i) are reduced to the single integrals 

Ac(t) = \ {J dT i a cosQ(n)^ , (51) 

A s {t) = l -(^J d n a (n) sinning . (52) 
This can be done since they have the structure 

A c ,s(t) = f dr f CyS (r)/ c , s (r) = i/ c 2 s (t), (53) 



where / CjS (r) is 

/c.(*) 



T ira(r)M, (54) 
7_ OG tsmfi(r)J 



and the index c stands for the cosine function and index s for sine. Integrating these integrals 
by parts, we obtain 

Z d (t) = 1 - - [(ai - a 3 + a 5 + ...) 2 + («2 - «4 + «e + --O 2 ] > 

Equations (46,47) and (55) give the probability amplitudes of the dark, bright and 
common states. They coincide with those one obtains if the successive transformations 
S n S n -i...SiS to the new set of d n b n c n -states are performed, as was done by Fleischhauer 
and coauthors in Ref. [37], i.e., |$ n ) = S n S„_i...SiS |$ >, where |$ ) = C x \1)+C 2 \2} + C 3 |3) 
is the initial state [C\{— oo) = 1, C 2 (— oo) = C^{— oo) = 0]. The transformation S is defined 
in Eq. (17) and the other transformations (Si) are specified below. These transformations 
have a simple meaning. Since the d — b and b — c transitions are excited simultaneously by 
the a and B "fields" , one can make a Si-transformation to the new set of dark, bright and 
common states, di, b\ and ci, where d\ is a particular mixture of the former d and c states, 
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h\ is a state orthogonal to state di, and C\ coincides with state b. The new mixing angle 
is aii = arc tan(d /B). Repeating this procedure n-times, one can get our solution if the 

conditions arc tan(d„ / B) ~a„ /B and \J a 2 n +B 2 pa B are applied at each step. 

Figure 4 (a-c) shows the comparison of the numerical solution (solid lines) with the 
expansions given by the equations (46), (47) and (55) (dashed lines). The parameters of the 
pulse train are specified above and they are the same as in Fig. 3 (a,b). Only the first two 
terms of the expansions are taken into account for each plot, which is justified because \a±\ >> 
\ot2\ ^> ^> \oti\ ^> .... For the X c and Y\> components, the maximum absolute values 
of the second terms of the expansions (i.e., the adiabatic terms) are already comparable 
with the amplitude of the oscillations (i.e., the nonadiabatic contribution, coming from 
the summation of an infinite number of the expansion terms). Of course, nonadiabatic 
oscillations are not described by the main part of the adiabatic solution presented by a few 
leading terms of the expansion. 

Concluding this section, we refer to a particular case when B(t) and a have the same 
time dependence. This is again the case of matched pulses (see the beginning of Section 
III), however in the dbc basis. Fleischhauer and coauthors, Ref. [37], classify this case as 
second order matched pulses. The solution of the equations (27)-(29) is trivial since these 
equations in terms of a new variable ( = F{r)dr can be reduced to a set of differential 
equations with constant coefficients, where F(t) = B(t)/B(to) —a (t)/ a (to). For the 
first time, this analytically solvable model was considered by Vitanov and Stenholm in Ref. 
[50]. We would classify the case as nonadiabatic, however in the second order d^ci-basis, 
where the transition takes place. If the generalized pulse area is properly chosen in this 
basis, the population transfer between the diabatic states 1 and 2 is complete. For these 
particular pulse areas, there are no nonadiabatic corrections, which is typical for the resonant 
nonadiabatic transitions. 

VI. NONADIABATIC CORRECTIONS 

All adiabatic terms tend to zero at t — > +oo, which secures for the three-level atom the 
adiabatic following of a particular state coinciding with state d at t — > +oo. However, as it 
will be shown below, if we sum all these infinitely small terms, the result will be finite. The 
net value of the small contributions of each adiabatic term is a nonadiabatic contribution, 
which specifies the excited probability amplitude left by the pulse train. To estimate this 
value, we rewrite the solution, Eqs. (42) and (43), as follows 

Y b (t) = f c {t) cosfi(t) + f s (t) sinfi(t), (56) 

X c (t) = f c (t) sinft(t) - f s (t) cosft(t). (57) 

When t — > +oo, the function f c (t) tends to a finite value whereas f s (t) tends to zero since 
in the corresponding integrals [see Eq. (54)] a (r) is an even function of time and f2(r) 
is an odd function. As a result, at the end of the pulse train Yj,(t) and X c (t) oscillate as 
cosfi(t) and sinf2(t), respectively, and they have constant amplitudes / c (+oo). The value 
of / c (+oo) defines the probability amplitude left by the pulse sequence in states b and c. 
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The excitation process is adiabatic if the nonadiabatic part / c (+oo) is small. If the 
nonadiabatic part becomes comparable with the main term ai(to) of the adiabatic expan- 
sion, then we can not describe the excitation process as adiabatic. Figure 5 (a) shows the 
numerically found dependences of the ai(t ) term (dashed line) and the nonadiabatic con- 
tribution / c (+oo) (solid line) on the maximum amplitude B of the pulses for the pulse train 
with Tr = 5. The adiabatic and nonadiabatic parts of the solution become comparable if 
B < 15r. 

The semilogarithmic plot of / c (+oo) versus B Q , Figure 5 (a), clearly demonstrates that 
the nonadiabatic part decreases exponentially with the increase of B . Laine and Stenholm 
[36] also found an exponential decrease of the nonadiabatic deviation from the ideal popula- 
tion transfer with increase of BqT (Bq is the amplitude of each pulse at the maximum and 
T is the interpulse distance, fixed in our case by the relation Tr = 5, so Bq/t is a variable). 
They followed the calculation procedure proposed by Davis and Pechukas [38] employing the 
Dykhne model [24]. According to Ref. [38], one has to find the zeros of the function B(t) in 
a complex plane t and take the one, t c , that is nearest to the real axis. Then the population 
of the dark state after the pulse train is 



|<i(+oo)| 2 oc exp 



-21m / B(t)dt 

Jt Q 



(58) 



We propose another calculation procedure of the nonadiabatic deviation. It will be shown 
that there are two nonadiabatic contributions, one coming from the Rabi frequency B(t) and 
another from the mixing parameter derivative a (t). Only the cooperative contribution of 
both determines the net nonadiabatic correction, while the Pechukas-Dykhne theory, taking 
into account only the _B(t)-change, underestimates the nonadiabatic correction. 

To show this, we express f c (t) via the Fourier transform of a (t): 

/+oo 
dra(T)e iuJ \ (59) 
-oo 

/t i r+oo 

dr cos Q(t)— / duo a{u) e~ iuJT . (60) 
■oo ^ J -oo 

Let us consider first the case if the Raman Rabi frequency is constant, i.e., B{t) = (3o = const. 
Then Q(t) = for and 

/•+OC l 

/ c (+oo)= / du-[6(u + p ) + 6(u-Po)]a(u) = a(u)\ u= i3 , (61) 

J-oo * 

where S(x) is the Dirac delta function and a(u) is an even function of uo. The Fourier 
transform of the mixing parameter derivative for a secant hyperbolic pulse train can be 
found, for example in Ref. [51]. This function is 

= sinh[far C taxL(smh(rT))] 
1 ^ sinh(^) 1 ^ 

If we take /3 = B(t ) = V2B sec h (rT/2), which is the value of the Raman Rabi frequency 
at time t — t when a (t) takes its maximum, then for our numerical example specified 
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above we obtain the amplitude of the nonadiabatic contribution / c (+oo) = 1.263 x 10~ 3 . 
This value is four times smaller than the amplitude of the nonadiabatic oscillations on the 
right tail of the functions X c (t) and Y b (t) (which is ~ 5 x 1CT 3 ), shown in Figures 4 (a,b), i.e., 
four times smaller than the nonadiabatic contribution given by the numerical calculations 
of the solution of Eqs. (27)-(29). 

To explain this difference and clarify the origin of the nonadiabatic contribution, we recall 
the interaction Hamiltonian in the o?6c-basis, Eq. (25) [see also Eqs. (15) and (26)]. This 
Hamiltonian resembles the interaction representation Hamiltonian of the three-level system 
excited by two resonant "fields" with amplitudes B(t) and a (t) (see Fig. 2). Assume, first, 
that the couplings B(t) and a (t) are absent and the system is in state d. Then, states 
d, b and c can be considered as having the same energies in the interaction representation. 
Switching on the coupling a (t) = a@(t) [here G(t) is the Heaviside step function and a is 
an arbitrary constant] mixes states d and b or in other words induces the transition d — > b. If 
the -B-field is also present and its amplitude (3 is constant, i.e., B(t) = p Q(t), this -B-field 
mixes states b and c producing a new couple of states b' and c', which are the states |l) and 

1 3) [see Eqs. (8) and (10)]. This couple is split by the energy gap 2(3 = 2B = 2^B~f+ B\ 
[see Eq. (6)]. This is the so called Mollow splitting [52] or quasi-energy splitting [46], 
[47]. Assume that without the -B-field the d-field is in resonance with levels d and b. The 
switching on of the .B-field mixes levels b and c, producing an additional splitting. Level 
b moves on the frequency B out of resonance with the a-field. If the a-field would have 
a 5-spectrum (in the case specified above it has only a zero frequency component), then 
it would not interact with the atom. However, because of the finite spectral width of the 
a-field, its spectrum has a component with frequency B on the far tail which is in resonance 
with the new position of level b. Only this spectral component excites the atom if the -B-field 
is on. With increase of the -B-field, the component of the a-field spectrum, which interacts 
with the atom, shifts further to the tail of the spectrum. If the .B-field amplitude changes in 
time, several spectral components of the d-field interact with the atom since at each instant 
of time some particular spectral component is in resonance. The process of the sweeping 
of the splitting B(t) along the tail of the d-field spectrum involves a broader band of the 
a-spectrum in the interaction. To find the net atom excitation in this case, we have to 
calculate the integral / c (+oo) where the change of B(t) is taken into account. This is done 
in the next two sections. 

Concluding this section, we apply our spectral method to find the approximate solution 
of the Rosen-Zener model. The comparison of our solution with the exact analytical solution 
helps to find the condition when our approximation is valid. Referring to the Bloch equations 
(30)-(32) discussed in Section III, we define the splitting of the two-level system as A = const 
and the coupling field amplitude as x — 2Vo sec h(rt), where Vq is the amplitude of the secant 
hyperbolic pulse of width r. According to the exact solution (see, for example, Ref. [10]), 
the population of the excited state e after the pulse is 



if initially the system is in the ground state g, i.e., p gg (—oo) = 1. 

Employing the correspondence between the Bloch equations for the two-level system and 
the Schrodinger equation for the three-level system discussed in Section III, we make the 




(63) 



17 



substitution a— \ an d B = A in Eq. (54) where Q(t) becomes a simple function of time, i.e., 
Q(t) = At. Then, if t — > +00, the solution is Yb(t) = f c (+oo) cos At, X c (t) = f c (+oo) sin At 
[see Eqs. (56)-(57) and the discussion following these equations]. / s (+oo) = since Q(t) is 
the odd function. / c (+oo) is simply the Fourier transform x(A) [defined as a nondimensional 
value according to Eqs. (59) and (61)] of the function x(^) at the frequency B = A, i.e., 
/ c (+oo) = x(A), [see Eqs. (61) and (59)], or explicitly 

/+00 
dt sec h(rt) cos At. (64) 
-00 

This integral also can be found in Ref. [51]: / c (+oo) = 27r(Vo/r) sec /i(7rA/2r). 

The population of the excited state e of the two-level system can be expressed as follows 
p ee = (l—w)/2. Since there is a reciprocity between the Bloch-vector components of the two- 
level system and the state probability amplitudes of the three-level system (w = Zj, v = 
u = X c ), the excited state population is p ee = (1 — Z c i)/2. The probability amplitudes are 
normalized, X 2 C + Y b 2 + Zj = 1, and hence Z d = a/ 1 — X? — Y b 2 = a/1 — / c 2 (+oo). Assuming 
that / c (+oo) < 1, we obtain 

p„( + oc).(^)\ectf(^). (65) 

Thus, our approximate result, Eq. (65), coincides with the exact solution, Eq. (63), if 
the population change is small and the sine function can be represented by the first term 
of its expansion in a power series. These conditions specify the validity of the adiabatic 
approximation. What is remarkable here is that the dependence of the excited state pop- 
ulation on the two-level splitting A in the approximate solution coincides with that in the 
exact solution. This shows that our spectral approach in the description of the nonadiabatic 
corrections is valid. 

If A — > 0, the nonadiabatic transition between states g and e takes place according to the 
straightforward solution of the B loch- equation for the two-level system: p ee (t) = sin 2 [6'(t)/2], 
where 6(t) = J_ x(j)dT. At t -> +00 we have p ee (+oo) = sin 2 (7rVo/r), which reproduces 
the exact Rosen-Zener solution (63) if A — > 0. 

VII. NONADIABATIC TRANSITION AS A QUANTUM JUMP: BASIC 

ARGUMENTS 

If the Raman Rabi frequency changes during the development of the mixing parameter 
a (t), those Rabi frequencies sweeping the frequency bandwidth of a (t) contribute to the 
nonadiabatic corrections. To take this process into account, we have to convolute the spectral 
content of both the mixing parameter and the Raman Rabi frequency. In general, this 
calculation is nontrivial. For the case of secant hyperbolic pulses [see the time development 
of the pulses, the Raman Rabi frequency and the mixing parameter in Fig. 3 (a,b)], it is 
possible to simplify the problem. Since B(t) has a minimum value of (3q at to where a (t) 
has its maximum, only the spectral components of a{uS) with ^ (3q contribute. Because 
time and frequency domains are inversely proportional, the main part of the nonadiabatic 
contribution appears in a short time interval around t , corresponding to the far tail of 
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the spectral distribution a(uo). This means that the time scale of the nonadiabatic change 
is shorter than while, according to the condition imposed on the parameters, the 

time variation of a (t) is much longer than the Rabi oscillation period ~ 1/A)- Thus, 
a nonadiabatic transition takes place almost jumpwise compared to the time scale of the 
variation of a (t). Due to this circumstance, we can simplify the calculation of the integral 
/ c (+oo) by expanding B(t) in a power series of t near t and retaining only the first two 
terms of the expansion 



B(t)^p [l + g(t-t ) 2 ], 



(66) 



where g = r 2 [3tanh 2 (rT/2) — l] /2. We verified this approximation by comparing numeri- 
cally two integrals / c (+oo) and /m(+°o) calculated with B(t) and its approximated value, 
Eq. (66), respectively. The comparison is shown in Fig. 5 b, where / c a(oo) (dashed line) 
is the approximation. The dependences of both integrals on the amplitude -Bo are indistin- 
guishable. 

The approximation of B(t) by a parabolic function helps to express / c (+oo) via the Airy 
integral. Further, to simplify the notations we set t = 0. Then the phase Q(t) is 



(67) 



and /ca(+oo) is expressed as 



fcA (+oo) = J dr cos (r + |r 3 )] J 
Evaluating the time integral in Eq. (68), we obtain 



dw a(u)e 



(68) 



/+oo 
dua(u)Ai [7(^0 + w)\ , 
-00 



(69) 



where Ai(x) is Airy integral and 7 = 1/ y/gfo. We derived this equation assuming that a(u) 
is an even function. The Airy integral has a different dependence for positive and negative 
arguments (see, for example Ref. [53]). If (3 + uo > 0, we have 



2 (Po+uy 

3 V 9Po 



2 m±uf 

3 V 9(3o 



(70) 



where I±±(x) is the modified Bessel function of order ±|. If the argument is negative, 
Po + uj < 0, then 



where J±i(x) is the Bessel function of order ±|. For large arguments both functions have 
a simple asymptotic behavior 
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cxp 



2 /WojW 

3 V flA> 



1 Ai [ 7 (/3 + w)] « — v (72) 

2v^[^o(/5o + ^)] 3 



if Pq + uj is positive, and 



cos 



2. /jffo+H 



3 



3V 9A 4 

7^ [ 7 (A, + w)] « — A (73) 

v^A^o + ^l) 3 

if Po + uj is negative. If the argument is zero, Po + uj = 0, this function is 

1 Ai [ 7 (/3 + u;)] U_ A = 3 r ^L- , (74) 

27rVv3#A) 

where T(|) « 2.679 is the gamma function. 

VIII. RABI CHIRPING ONLY 

If the a-field does not depend on time and has a value a (t) —a (to), then a(uj) = 2n a 
(to)8(uj)/r and the main contribution to the nonadiabatic part is given by the Airy integral 

/ cA (+oo) = 2Tt 1 Ai{ 1 p Q ) a (t )/r, (75) 

which has the explicit form 

f cA (+oo) = tanh ,/^-Iexp ( , (76) 



expressed via approximation (72). In this section we show how the approximation, specified 
above, is related to the Dykhne-Pechukas model [24], [38]. 

For large Po, the Airy integral can be calculated by the saddle point methods as shown, 
for example, in Ref. [54]. First, the stationary or saddle point t s is found where the phase 
Q(t) becomes stationary: Q (t s ) = 0. In this point, the Raman Rabi frequency becomes zero 
since Q (t) = B(t). For the case of the positive argument of the Airy integral, the saddle 
point is in the complex plane and t s has only an imaginary component [Re(t s ) = 0], i.e., 



ts r]] 3tanh 2 (rT/2) - l' ^ 

In fact, there are two saddle points: t s = rki/^fg. The negative sign is chosen to avoid 
exponentially increasing numbers. Then the method of stationary phase (one of the saddle 
point methods) is applied. This method is applicable to integrals of the form 

I(P ) = [ e lM{z) dz, (78) 
Jc 
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where (3q is large and C is a path in the complex plane such that the ends of the path do not 
contribute significantly to the integral. The idea of the method is to deform the contour C 
so that the region of most of the contribution to I(/3o) is compressed into as short a space 
as possible. This compression occurs at the saddle point. 

Since the main contribution to this integral comes from the vicinity of the saddle point, 
one may conclude that in the adiabatic limit the atom abruptly makes a nonadiabatic 
transition between adiabatic states. This is very much like crossing a Stokes line for the 
asymptotic series, and resembles the turning point problem in the WKB theory (see, for 
example, Refs. [10], [54]). In the WKB method, the one-dimensional Schrodinger equation 
is considered for a particle with mass m in a potential V(x): 

d 2 ^ 2m rrn Tr/ . -. 

_ = __[£ _„(,)]». (79) 

The turning point is the one where the energy of the particle coincides with the potential: 
E — V(x) = 0. If 2m[E — V(x)]/h 2 oc x, the solution of the equation (79) can be expressed 
via the Airy integral [54]. In this approach the WKB connection formula, relating the 
exponentially small solution on one side of a turning point to an oscillatory solution on the 
other side, is derived quite naturally. 

The Pechukas-Dykhne [38] recipe of the calculation of the nonadiabatic contribution 
is similar to the method described above since the definitions of the saddle point t s and 
the crossing point t c , where B(t c ) is zero, are the same. However, in our consideration we 
substituted the expression for the Raman Rabi frequency 

= 2 Vl + co,h( r r)cosh(2rt) 
V ; cosh(rT) + cosh(2rt) ' V ; 

by the expansion near time to, Eq. (66). Therefore, we have only two saddle or crossing 
points, (79), whereas equation (80) has an infinite number of crossing points in the complex 
plane. According to Pechukas, one has to take the crossing point that is nearest to the real 
axis. For the secant hyperbolic pulse train, this point, as shown by Stenholm [36], has only 
an imaginary part lm(t c ), i.e., 



i 

t c — - tan 
r 



coth 



tT 
~2~ 



(81) 



Reformulating the Dykhne approach [24], we conclude: if there are no spectral components 
of a(u) matching the frequency gap between the quasi-energy levels split by the 5-field, the 
d-field comes to resonance with the d — c transition at the imaginary time t c when B is zero. 

This approach disregards the spectral content of a (t). If we take a (t) =a (t ) through- 
out the excitation, this method gives also an underestimated value of the nonadiabatic con- 
tribution. For our numerical example, equation (76) gives / c a(+°o) = 7.131 x 10~ 4 , which 
is 7 times smaller than the value of the nonadiabatic contribution given by the numerical 
calculation of the Schrodinger equation (27)-(29). Figure 6 (a) shows the comparison of the 
numerically calculated nonadiabatic component / c (+oo) (solid line), Eq. (54), with that 
calculated for the case if the time dependence of the a (t)-field is neglected (dashed line). 
The latter is given by Eq.(76). The plots obviously demonstrate that this approximation 
underestimates the nonadiabatic contribution. 
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Concluding this Section, we demonstrate the application of our method to the Landau- 
Zener model. In this model the parameters of the corresponding two-level system are A = 
2Xt and \ — 2V . Here, the level crossing takes place ait — and the nonadiabatic correcton 
appears in the vicinity of the crossing point. After the substitution specified in Section III 
and repeating the same procedure as described at the end of Section IV, we find that 

/+oo 
dtcos(Xt 2 ), (82) 
-oo 



/+oo 
dtsin(Xt 2 ). (83) 
-oo 

The value / s (+oo) is not zero for this model since Q(t) = Xt 2 is even function of time. The 
integrals in Eqs. (82)- (83) are expressed via Fresnel integral asymptotes (see, for example, 
[53]): / c (+oo) = / s (+oo) = Vo^F/X. Recall that Z d (+oo) = y/1 - f 2 (+oo) - / s 2 (+oo). 
Then in the adiabatic limit, V 2 <C X/2tc, we obtain 

. l-Z d (+oo) Vq 
pU+oo) = | ^ tv^-. (84) 

This expression coincides with the leading term of the exact solution, Ref. [10], expansion, 
that is 

/ T/2\ T/2 

p ee (+oo) = 1 - exp f -7T^ J ~ 7T^. (85) 

We have to emphasize that our solution is valid only if the atom adiabatically follows 
the ground state g. Therefore, the nonadiabatic correction (84), specifying the population 
probability of the excited state e, coincides with the Landau-Zener solution, Eq. (85), if 
Vq <C X/2n, i.e., when the population of the e and g states does not change appreciably. 
In this case our three-level system dbc, equivalent to the corresponding two-level system 
eg, also follows adiabatically state d. The opposite case is considered in the Appendix. 
There the new adiabatic states d\b\C\ are introduced whose population also does not change 
appreciably. This transformation is necessary if Vq ^> A/27T, because the reciprocal three- 
level system dbc does not follow state d, making a transition from state d to state c and back 
at the level crossing. 

The necessity of the transformation can also be formulated in a different way. To apply 
the adiabatic solution, we need a certain relation between the interaction parameters B and 
a. Then the expansion in a's, where each next term of the expansion is smaller than the 
previous one [see Eqs. (42)-(55)], is applicable. If V 2 <C X/2-n, this expansion can be made 
in the <i&c-basis. In the opposite case, this expansion is not applicable. Therefore, we have 
to find new dibiCi-b&sis where the new interaction parameters B and a x are related such 
that the expansion in a's converges. 
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IX. COOPERATIVE CONTRIBUTION OF RABI CHIRPING AND THE TIME 
DEPENDENT COUPLING, THE APPLICATION TO STIRAP BY SECANT 

HYPERBOLIC PULSES 



In previous sections we showed that the contribution to the nonadiabatic corrections 
of the Rabi frequency chirping or the time dependence of the a (t)-field, being taken into 
account separately, was not enough to have the right nonadiabatic correction for STIRAP 
induced by secant hyperbolic pulses. Only both processes taken together, i.e., the Rabi 
frequency B(t) chirping and the participation of all spectral components a(u>) of the a (t)- 
field engaged by this chirping to excite the atom, give the right value of the nonadiabatic 
component. This value is given by Eq. (69). To calculate analytically the convolution 
integral of the two spectra, a(u) and the Airy integral describing the process of the Rabi 
frequency sweeping, we make two approximations. First, we take the approximation of the 
Airy integral, given by Eq. (72) in the frequency domain u > — /3 - This gives a small 
overestimation of the integrand near — (3q. For this reason, we start the integration from this 
value, not from — oo. The oscillating part of the Airy integral for oo < —(3q gives a much 
smaller contribution than the main part, which is located between — (3$ and zero, (— /3o,0). 
So, neglecting the part (— oo, — /3q), we compensate the overestimation near — /3 . Second, 
we approximate the spectrum a{uS) in the domain (— /3 ,0) by 



where R = (n — arc tan[sinh(rT)]) /2r ps n/Ar. This approximation also gives a slight 
overestimation of the integrand near u ~ 0. To compensate this, we stop the integration 
at u — 0. The integrand has a maximum between u = —/3q and u — 0. To calculate the 
contribution of this part, we use the modified method of the saddle point method. Usually, 
in the method of the saddle point, the integration near the point is extended to ±oo. In 
our case, to avoid the overestimation of the integrand we limit the integration by finite 
boundaries. The calculation of these boundaries for the deformed integration contour C is 
simple because the deformed C stays on the real axis. The result of the integration is 



where Q is a correction factor, which takes into account the finite integration boundaries to 
avoid the overestimation of the integral. Explicitly, we have 



a{uj) = n exp (-R , 



(86) 




(87) 



Q = - [erf(/i max ) + erf(/i min )] , 



(88) 




(89) 




(90) 
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and erf (x) is the error function. Figure 6 (b) shows the comparison of the true nonadiabatic 
contribution / c (+oo) (solid line) with our approximate calculation of / c a(+oo), given by Eq. 
(87) (dashed line). Figure 6 (c) compares the same dependencies if the correction factor is 
Q — 1 (dash dotted line). 

For our numerical example specified above, the approximated value of the nonadiabatic 
contribution is f cA (+oo) = 5.822 x 10~ 3 . As was discussed above, this contribution appears 
in a very short time interval around to- Therefore, we can approximate the solution of the 
Schrodinger equations (27)-(29) by 

X bA (t) ^a 2 -a 4 + e(t - t ) f cA (+oo) cosft(f), (91) 

X cA (t) ^ ai -a 3 + 0(f - t ) / c a(+oo) sinft(t). (92) 

where the index A stands for the approximation, 0(t — to) is the Heaviside step function, 
Q(t) = Jq B(t)cIt and B has its exact value [not the approximation (66)]. Figures 7 (a,b) 
show the comparison of the numerically found solutions of the Schrodinger equation (solid 
lines) with the approximation (91)-(92) (dashed lines). The fit of the solutions is striking. 
This means that the nonadiabatic contribution really appears in a quite short time range 
around to- 

X. CONCLUSION 

The introduction of the basis of the dark and bright states facilitates the understanding of 
the physical processes in the three-level atom excited by a bichromatic field. The dynamic 
evolution of the atom between the dark, bright and common states is as simple as the 
dynamics of the two- level atom excited by one field with Rabi frequency B. This is because 
the evolution of the three-level atom can be effectively reduced to the evolution between two 
states, i.e., the bright and common states. The reduction of the three-level model to the 
two-level one allows also the application of the Bloch- vector model and Bloch equations for 
the treatment of the three-level atom excitation by the bichromatic field. 

This quite simple algebra is applicable for the case of matched pulses. If the pulses do not 
match in shape and have different time dependencies, one can also reduce the consideration 
to the Bloch-vector model since there is a similarity between the Schrodinger equations for 
the probability amplitudes of the dark, bright and common states and the Bloch equation for 
an effective two-level system. The effective detuning of the two-level system from resonance 
is B(t) and the Rabi frequency is a (t), which is the derivative of the mixing angle in the 
dark state development in states of the reciprocal three-level system. This similarity allows a 
simple interpretation of the physical processes in the three-level system in case of adiabatic 
following of the dark state. In terms of the Bloch-vector model, the adiabatic following 
criterion is formulated resulting in three conditions. First, the initial state must coincide 
with the Bloch-vector aligned along an effective field at t = — oo. This field is determined 

by the value \J B 2 (t)+ a 2 (t). Second, the effective field makes a small angle with its initial 
position when the reciprocal parameter a (t) responsible for the atom-field interaction takes 
its maximum. Third, the time scale of the change of the "pulse", determined by the scale 
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of the variation a (t), is much longer that the period of oscillations induced by the effective 
field B(t). If all three conditions are met, the dynamic evolution of the three-level atom 
excited by a sequence of pulses is described by simple algebra. The adiabatic population 
transfer has a nice interpretation in terms of dark, bright and common states. 

We developed an approximation describing the adiabatic interaction of the three-level 
atom with two resonant pulses. The method of estimation of the nonadiabatic correction 
is presented. It is applied to the case of two secant hyperbolic pulses. The adiabatic part 
of the solution describes the excitation and de-excitation processes of the three-level atom. 
The time dependence of the adiabatic part is smooth and follows the derivatives of the 
mixing angle a(t). Both parts, excitation and de-excitation, are symmetric in time with 
respect to t — where a (t) takes its maximum and they exactly compensate each other. In 
this respect the adiabatic following resembles the soliton-like interaction with the field. The 
nonadiabatic part appears in a short time interval in the vicinity of the maximum of the 
mixing parameter a (t). It contains the information about the excitation left in the atom 
by the pulses because of the imperfect following of the dark state. Since this excitation 
lasts only a short time, we consider the nonadiabatic process as a transition between the 
ground and excited states, which takes place like a jump. The origin of the transition has a 
simple interpretation. The parameter B defines the coupling strength of the bright (b) and 
common (c) states of the three-level system in the ci&c-basis. The parameter a (t) defines the 
coupling strength of the bright and dark (d) states. It is assumed that, initially, the system 
is in the dark state and any time we have B ^>a (t). The 5-coupling moves the bright 
state from the resonance with the a (t)-coupling to the value determined by the frequency 
B. If the spectral content of the a (t)-coupling has a component with frequency B, the 
nonadiabatic transition takes place. If B changes in time, several spectral components of 
the a (t)-coupling contribute to the transition. Our results are compared with those of 
Laine-Stenholm [36] and Fleischhauer et al. [37]. 

We applied also our method to the Rosen-Zener and Landau-Zener models. The com- 
parison of our approximate solution with the exact solutions available for these models is 
quite encouraging. We obtained exponentially small nonadiabatic corrections for the case of 
a symmetric time dependence of the excitation parameters B(t) and a (t) with respect to 
t — (our model, the Rosen-Zener model and the Landau-Zener model with Vq ^> \/2ir), 
and a small power dependent nonadiabatic correction for the case of an asymmetric time de- 
pendence of B(t) (Landau-Zener model with V 2 <C \/2ir). We considered also the excitation 
of the three-level system by asymmetric pulses and found a power dependent nonadiabatic 
correction. The results on asymmetric pulses, exciting the three-level atom, will be presented 
in another paper. 

The simplified algebra developed in this paper could be useful for the description of the 
atom state manipulation by coherent fields. 
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XII. APPENDIX 



In this section we compare our approximate solution with an exact solution of the 
Landau-Zener model if the parameters of the corresponding two-level system A = 2 At and 
X = 2V , satisfy the condition V 2 3> A/27T. This condition is opposite to the one considered 
at the end of Section VII. There, since the coupling Vq is small compared to the sweeping 
rate A of the two-level system splitting A (V 2 <C \/2ir), the system adiabatically follows 
the ground state g and the net population of the excited state e is small for t — > +00. If 
Vq ^> A/27T, states g and e are not an adiabatic basis. To apply the adiabatic approximation 
in this case, we have to choose a new basis that is adiabatic. This can be done if we use 
again the similarity between the two-level and three-level systems discussed in Section II. We 
consider the three-level system dbc described by the Hamiltonian (25) with the parameters 
B = A = 2\t and a— x = 2V . Then we make a Si-transformation to the new basis dib\C\ 

\di) = sin cci \d) + i cos a± |c) , (93) 
\bi) = i cos ai \d) + sin «i |c) , (94) 
|ci> = \b) , (95) 



where a\ = arc tan (A/x) and the new Rabi frequency is B — ^/ A 2 + x 2 - A new mixing 
angle a\ takes the values ±7r/2 at t = ±00. Thus, it varies from — tc/2 to +ir/2 with time 
and, hence, state d\ coincides with — \d) and + \d) before and after the crossing of levels g 
and e, which takes place at t — 0. If the system adiabatically follows the changing state d±, 
no transition to the b\ state takes place and, finally, the system is left in the initial state d 
but with the 7r-shift of its phase. This phase shift originates from the transition d — > c — > d 
in the ci&c-basis. One can find that, according to Eq.(93), the system is in state d at t = ±00 
when a\ = ±7r/2, and in state c at t = when ot\ — if the system follows state d±. 

We know from the exact solution of the Landau-Zener model that p gg (+oo) = 
exp(— tiVq/X) and p ee (+oo) = 1 — exp(— hVq/X) if p gg (—oo) = 1 [see Eq. (85)]. Thus, 
if ttVq 3> A, the population of the ground state g of the corresponding two-level system is 
exponentially small at t — > +00. With the help of our method for the estimation of the 
nonadiabatic correction, we can find the probability amplitudes of states b\ and c\ left by 
the interaction x (coupling states d and b) after the level crossing of the equivalent two-level 
system. These probability amplitudes allow to estimate the extent of the depopulation of 
state d\. If the system adiabatically follows the d\ state, the state- vector changes from — \d) 
to + \ d). Since the development coefficient Cd of the state vector in the <i&c-basis corresponds 
to the Z-component of the Bloch-vector of the equivalent two-level system, the change of 
the Cd sign means the flip of the Z-component. This flip means a population inversion of 
the two-level system, i.e., if before level crossing (at t = —00) the system is in the ground 
state g, after the level crossing (at t = +00) it is in the excited state e. We can estimate 
the amount of population left in the ground state, which is the nonadiabatic correction. 

To find the nonadiabatic correction, we apply the 5i-transformation 



Si = 



sinai — icosai 
-icosai sinai 




(96) 
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to the Hamiltonian (25) with the parameters B = 2Xt and a= 2V . The result is 



H-diha — Hdibia + i Si S 1 1 , (97) 

where 

n dlblC1 = -V^ 2 + x 2 (n ici + P Clbl ) , (98) 

and 

% Si S± l = - Oi X (P bldl + P dl6l ) . (99) 

Here, for simplicity, we set the phases [see Eq. (1),(4)] tpi and </?2 equal to zero. The 
development coefficients C dl , C bl , C C1 of the state vector \&dibia) i n the (i^Ci-basis satisfy 
equations that are similar to equations (27)-(29) if the substitution Zi = — C dl , Y 1 = iC bl 
and X 1 = C C1 is made. These equations are 

Zi= - Oi X Y u (100) 

Yi= -BX 1 + dx Zi, (101) 

Xi= SFi. (102) 
The derivative of the mixing parameter a\ has an explicit form 



A 1 

Vn 



«i= v 7-^2- ( 103 ) 



This function has a bell shape with maximum value X/Vq & t t = 0. The modified Rabi 

frequency B = 2Vq\J\ + (At/Vo) 2 has a minimum value 2Vo at t — and increases to 
infinity at t — ±oo. Qualitatively, the time dependence of the interaction parameters 
resembles what we have for the STIRAP induced by the secant hyperbolic pulses, i.e., B(t) 
and a (t). Therefore, we can apply the same procedure of the calculation of the nonadiabatic 
correction and the same approximations. 

First, we approximate the Rabi frequency by its expansion in a power series of (At/Vo) 2 , 
retaining only the first two terms 

S(t)«A>(l + flfc 2 ), (104) 

where (3q = 2Vo and g = (A/Vo) 2 /2 [see Eq. (66)]. Rescaling time as r = Xt/Vo, we can 
define the nonadiabatic correction / Cj 4(+oo) as follows 

f + °° cosfi(r) , , N 

f cA (+oo) = / \ J dr, 105 

where 



27 



(106) 



and T) = Vq/\, g± — 1/2 [compare the function fi(r) with the one in Eq. (67)]. The Fourier 
transform 



/+oo 
dr «i (r)e~ 
-oo 



of the coupling parameter «i (r) is 

ai(u) = 7rexp (— |cj|) . 
Following the same procedure as the one described in Section VII, we obtain 

/ + °° r j i 

ai(u)-q~^Ai rj~3 (2rj + u)\ du. 
■CO 



(107) 



(108) 



(109) 



Employing the same approximation of the Airy function by equation (72) for positive argu- 
ments and taking into account only the negative wing of the function ai(u), we limit the 
integration in Eq. (109) by the boundaries (— 2rj, 0). These boundaries are taken to reduce 
the overestimation of the integrand introduced by the approximation applied to the Airy 
function and ai(u) function. The arguments justifying this procedure are absolutely similar 
to those that are applied in Section IX to the case of the secant hyperbolic pulses. 
After some algebra, we obtain 



/r.i( 5c) ttQj <^p ( ~r) ) . 



(110) 



where 



and 



Qi = - [erf(/i max ) + erf(/i min )] , 



111^ 



\ 



2r] 



2V2 



2\ 3 1 
3 +6 



(112) 



(113) 



This result can be obtained from Eqs. (87)-(90) by the straightforward substitution R — 1, 
Po = 2r) and g = g 1 . 

We take 00) = 1 or C^— 00) = —1 and Cd(— 00) = 1 as initial condition, which 
means that the reciprocal two-level system is initially in the ground state g. The final state 
is defined by the relation Cd(+oo) = C^+oo) = — Zi(+oo). Since Zi keeps its sign, 
changes sign. The population difference of the equivalent two-level system w = p gg — p ee is 
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related to Cd as w = Cd- If Cd becomes negative, then p ee > p gg . In this case, to estimate 
the population fraction left in the ground state g, we have to use the expression 



l + C d (+oo) l-Zi(+oo) 



(114) 



After the substitution Z^+oo) = a/1 — /^(+oo) into Eq. (114) and expanding the square 
root in a power series, we find the nonadiabatic correction for the Rosen-Zener model 



p gg (+oo) w -/^(+oo) = 



(115) 



Figure 8 shows the comparison of our approximate solution, Eq. (115), with the exact 
solution of the Landau-Zener model 



(116) 



The difference is small and it originates from the approximations made. To estimate the 
accuracy of our result, we compare Eqs. (115) and (116) in logarithmic scale as follows 



In [p ee (+oo)] = -i]7i 



-V 



10 2 , 
In 

3 T] 



(117) 



For large r], the accuracy of the exponential factor calculation, expressed in the square 
brackets in Eq. (117), is 6%. 
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XIII. FIGURE CAPTIONS 



Fig.l. The excitation scheme of the three- level atom by two resonant fields with inter- 
action constants B x and B 2 : A-scheme (a), and the same excitation scheme in the basis of 
dark d, bright b and common c states (b). B — \J B\ + B 2 is the generalized interaction 
constant for the Raman excitation. The vertical scale in diagram b is not the energy and 
represents the initial population of the levels before the excitation. The higher populated 
level is put at the bottom and the lower or unpopulated level is at the top. 

Fig. 2. The bichromatic excitation scheme of the three-level atom in the dbc basis if the 
state developments of \d) and \b) change in time. The coupling parameter of states d and b 
is a, the derivative of the mixing angle (see the text). 

Fig. 3. (a) The pulse train with corresponding interaction parameters Bi(t) and B 2 (t) 
evolving in time (bold lines). They are normalized by the maximum value B Q . The delay 
between pulses is rr = 5, t = 0. The time is scaled in units of r. The time dependence of 
the mixing parameter a is shown by the thin line. The dashed line shows the dependence 
of the mixing parameter derivative a normalized by r. In (b), the time dependence of the 
bichromatic Rabi frequency %/2 = B(t) (solid line) and the mixing parameter derivative 
(dashed line) are shown for comparison. 

Fig. 4. (a) The evolution of the amplitudes X c (plot a), Y b (plot b), Z d (plot c) for a pulse 
train Bi(t), B 2 (t) with the parameters B = 42. 8r, t\ = 2.5/r, t 2 = — 2.5/r. Solid lines are 
the numerical solution of equations (27)-(29). Dashed lines are analytical approximations 
given by the first two terms of equations (46,47) and the first two terms in each of the 
parentheses of equation (55). 

Fig. 5. (a) Comparison of the dependence of the adiabatic ai(to) (dashed line) and 
nonadiabatic / c (+oo) (solid line) parts of the analytical solution of Eqs. (27)-(29) on the 
maximum pulse amplitudes -Bo- In (b), the comparison of the nonadiabatic part, numerically 
calculated using the actual Rabi frequency B(t) (solid line) and its parabolic approximation, 
Eq. (66) (dashed line), are shown. / c (+oo) is the true nonadiabatic part and / c a(+oo) is 
an approximation. 

Fig. 6. (a) Comparison of the dependencies of the nonadiabatic contributions / c (+oo) 
(solid line) and / c a(+oo) (dashed line) versus the pulse amplitude B . The first is calculated 
with and the second without taking into account the time dependence of the mixing param- 
eter derivative a (t). (b) The plots of the true, / c (+oo), and the approximate, / c a(+oo), 
nonadiabatic contributions versus B , where / c a(+oo) (long dashed line) is calculated taking 
into account the time dependence of a (t) using the approximation described in the text, 
(c) The same plots as in (b), except f C A(+oo) (dot dashed line), where the correction factor 
Q is dropped (see the text). 

Fig. 7. Numerical solution of the Schrodinger equation (solid lines) for the amplitudes 
X c (t) (plot a), Yb(t) (plot b) and the approximation given by equations (87)-(88) (dashed 
lines) for the amplitudes X cA (t) (plot a), Y bA (t) (plot b). 

Fig. 8. Comparison of the exact solution of the Landau-Zener model (dashed line) with 
our approximate solution (solid line) for the ground state population of the reciprocal two- 
level system after level crossing. The dependence of the exact solution of the Landau-Zener 
model for the ground state population, [p gg (r))]Lz, on the parameter rj = Vq /A is shown by 
the dashed line [see Eq. (116)]. The same dependence of our approximate solution, [p gg (r])]A, 
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is shown by the solid line. Our solution is given by Eq. (115). 
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